A SOLUTION TO AN OPEN PROBLEM ON THE MATHIEU SERIES 
POSED BY HOORFAR AND QI 



GERGO NEMES 



Abstract. In this paper, the author gives a solution to an open problem about the famous 
Mathieu series, that is, he obtains a sharp double inequality for bounding this series. 



1. Introduction 



We consider the series 



(1-1) S(r) = J2 -, '"2 , r> ° 

„>i (n 2 + r 2 ) 

called Mathieu's series in the literature. It was introduced by Mathieu in his 1890 work on 
elasticity of solid bodies. Several interesting problems, solutions and bounds for the Mathieu 
series can be found in the works we collected in the references. Alzer, Brenner and Ruehr [1] 
showed that the best constants a and j3 for which 

-5- — < S (r) < 1 - 

r 2 + a w r 2 + f3 

holds, are a — 21^(3) an< ^ & = j|> wnere C denotes Riemann's Zeta function. Our main purpose 
is to prove the following problem proposed by A. Hoorfar and F. Qi [5]: find the best possible 
constants a and b such that 

(1.2) J- <S(r)< =- 

9 1 4r 2 + 1 ! I 4r 2 + 1 

r 2 + _ __( r 2 +a ) r 2 + _ __( r 2 +fe ) 

holds true for all r > 0. They showed that a < 3/2 and 6 > 1/4. The following theorem gives 
the complete answer. 

Theorem 1.1. For any real number r > 0, the inequality (ll.2[) ZioZds tuif/i i/ie 6esi possible 

C (3) 13 

constants a = = 0.9915168156 . . . and b = — . 

6(, (3) — 6 30 

2. Lemmas 

To prove our theorem we define a function a for any real number r > by the equality 

1 



S(r) = 



yr £ + a (r)) 



2 12 

Theorem 11.11 is the corollary of the following lemmas 
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Lemma 2.1. For any real number r > 2.57, a(r) is a strictly decreasing function of r, satisfies 
the inequalities || < a(r) < g^$^~g ? and limr^+oo a (r) = 

Lemma 2.2. for any rea/ number < r < 0.3, a(r) is a strictly decreasing function of r and 

13 < QfM< c(3) 
30 "U/l 6C(3)-6' 



satisfies the inequalities < a (r) < — 



Lemma 2.3. For any rea/ number 0.3 < r < 2.57, a(r) satisfies the inequalities || < a (r) < 
CO) 



6C(3)-6- 

In this paper we use the following definition of the Bernoulli numbers: 

— - — = 1 - -x + B 2 ^-r - B^— + B 6 ^— , Ixl < 2-7T, x + 0, 

e x - 1 2 2! 4! 6! 11 ^ 

thus the first few are 

B2 = l> Be = 42' Bs= 3^' Sl0 = i' ■■■■ 

To prove our lemmas we shall use the following theorem of Russell [T^] . 

Lemma 2.4. Let S (r) be defined by (TJ\) and f (x) = xj (e x - 1) if x ± 0, / (0) = 1. Then 



5 W = i-EJ^+^W (r>0, fc=l,2,...), 
»=i 

where Rk (r) may 6e expressed in either of the forms 

i? fc (r) = ±> / ( x ) cos (ra) dx 



Nfe + l 







= -%rj o f (2k+2 Hx) S iMrx)dx 

and 

3. Proofs of the lemmas 
Proof of Lemma \2.1\ From the definition of a (r) we have 

(l2r 4 + 2r 2 -l)S(r)-12r 2 

(3 ' 1} " W = 12-(12r 2 + 6)S(r) ' 

According to Hoorfar and Qi we find that 

1 , N 3 
- < a(r) < -, 
4 " w - 2' 

if r > 0. We have to improve this for r > 2.57. Note that if 
(3.2) A(r) < S(r) < B(r) 

and 

(3-3) 2 ^ <A{r), B(r), 

r z + 1/2 

then 

(I2r 4 + 2r 2 - l) B (r) - 12r 2 (l2r 4 + 2r 2 — l) A (r) - 12? 



12 - (12r 2 + 6) B (r) w 12 - (12r 2 + 6) A (r) 
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From Lemma l2.4l we have 



1 



6r 4 30r 6 42r 8 



ir{7/2\ 1 
42 



1 1 



1 



1 



53 



> 



if r > 1. Let 



then 



(3.5) a(r)< 



6r 4 30r 6 500r 8 j- 
53 



1/2' 



6r 4 30r 6 500r 8 ' 



(l2r 4 + 2r 2 - l) G (r) - 12r 2 _ 2600r 6 + 1758r 4 + 268r 2 - 159 



< 



C(3) 



12 -(12r 2 + 6) G(r) 6000r 6 - 2100r 4 - 2208r 2 - 954 6C(3)-6' 

if r > 1.3. Applying Russell's theorem again, we obtain 



1 



1/2 



<*«<i- 1 



i 



i 



if r > 0. Let 



6r 4 30r 6 42r 8 



ifM = 1 L 



tt/7/2\J_ 
2 V 3 J 42 



1 1 



1 



1 



r 2 6r 4 30r 6 50r 8 ' 



then 



a (r) > 



6r 4 30r 6 50r 8 ' 
(l2r 4 + 2r 2 - l) K (r) - 12r 2 _ 260r 6 - 123r 4 - 23r 2 + 9 



12 - (12r 2 + 6) K (r) 
From this and from (13.51) we hnd 



lim 



260r 6 - 123r 4 - 23r 2 + 9 
*oo 600r 6 - 21 Or 4 + 78r 2 + 54 



< lim a (r) < lim 



600r 6 -210r 4 + 78r 2 + 54' 

2600r 6 + 1758r 4 + 268r 2 



159 



Since 



lim 



260r 6 - 123r 4 - 23r 2 + 9 
r^Too 600r 6 - 210r 4 + 78r 2 + 54 
we obtain 



lim 



r^+oo 6000r 6 - 2100r 4 - 2208r 2 - 954' 

2600r 6 + 1758r 4 + 268r 2 - 159 _ 13 
6000r 6 - 2100r 4 - 2208r 2 - 954 ~~ 30 



lim a (r) 



13 
30' 



Now we consider the monotonicity of a. We will show that the value of the derivative a' (r) 
is negative for any real number r > 2.57. A straightforward calculation gives 



2r z 



2 n 



(r 2 + n 2 y 



Thus it is sufficient to show that for r > 2.57, 

(3.7) T (r) = (V 4 + 3r 2 + ±) S 2 (r) - (6r 2 + 2) S (r) - (V + ^) E 



2n 



(r 2 + ?7 2 ) 3 



3 > 0. 



From Lemma l2.4l we have 
S(r) 

S 2 (r) 



1 



1 



1 



1 



6r 4 30r 6 42r 8 30r 10 



+ Ra (r) , 



2R 4 (r) Rj (r) - 3 R 4 (r) + 5 60R 4 (r) + 49 
~~F 2 3T 4 15r 6 1260r 8 

42i? 4 (r) + 23 



630r 10 



Af (r) 
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y 2n 1 , n = l J_ _L _L_ _^ iff'M 

l^( r2+n 2f 4r 1 J 2r 4 6r 6 20r 8 21r 10 12r 12 4r 4 ^ ' 



n> 

where 



121 4 41 11 2 

(3.8) M (r) - - 2100rl2 + + i4700r 16 + 630^8 + 900^° + 4 W ' 

By substituting these into (13. 7[) we obtain 

W W W 2 8r 5r 2 210r 4 210r 6 

282i? 4 (r) - 301 6i? 4 (r) - 43 1 6i? 4 (r) + M (r) 
1260r 8 180r 10 + 24r 12 + 2 ' 

From Lemma EU we have 



_ . . 5 691 5 691 1 7T /l3/2\ 



5 J_ 

66r 12 r 14 



Thus if r > 1, 



i? 4 (r) 


107i? 4 (r; 


5r 2 


210r 4 


43 


43 


180r 8 


+ 180r 10 + 


43 


43 


180r 8 


+ 180r 10 + 



1451 



> 
> 

This means that 

T (r) > 3r 4 M (r) + 3r 2 M (r) + 



210r 6 1260r 8 180r 10 

1 2237 6667 4979 3067 



, ,_2,,^ , ^(^ , ^4 (r) , 43 43 



2 8r 180r 8 180r 10 

1 1 1451 6i? 4 (r) + M (r) 



24r 12 66r 14 1155r 16 2 

From (13.81) it is clear that 

121 | 4 | 41 

W > 2100r 12 + 315r 14 + 14700r 16 ' 

We also have 

5 691 5 691 J_^A 3 / 2 A R 

4(r) " 66ri 2 2730r 14 + 6 (rj > 66r-i 2 2730r-i 4 r 44 2 ^ 6 1 

5 3967 

> 



66r 12 2730r 14 ' 
and thus 

i? 4 (r)r i^(r) 104 164 13579 27302953 405829 



(3 9) T (r) > 

v ' ' w 2 ' 8r ' 1575r 8 ' 1575r 10 80850r 12 6306300r 14 323400r 16 ' 
Now we are going to estimate i? 4 (r). We have 

^W=(-66^- 2730^ +i?6(r) =ll^ + 195^ + ^ (r) ' 
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where 



R' 6 (r) = 



/ (13) (x) cos (rx) dx 
14 f 00 1 f 00 

— — I / (13) (a;) cos (rx) dx — / x/ (13) (x) sin (rx) dx 

r Lb Jq r 14 7 

x/( 13 ) (x) sin (rx) dx. 



14 „ , N 1 
—^e (r) - -j 



Russell showed that 



{k) ■ fc!2-C(fc) 



f {k) (27TX) 



27! 



(2 ^^ i(n2 + , 2) (^)/ 2 ' 

if fc > 2. Diananda [3] proved the inequality 

2n 1 



E 



< 



Thus 



' (n 2 +x 2 f l)x 2 ^- 2 
fc! 2 



, /i > 1, x > 0. 



/ (fc) (27TX) 



From these we obtain 



{2ir) k (k-l)x k 

\R«(r)\<—\R a (r)\ + ± x /( 13 )(x) |sin(rx)|dx 
r r 14 J 



j, k > 1, x > 0. 



< 



< 



4tt 2 r 00 


x 


14 1 7T 


(T) 


2 




1051 


49 


100r 14 H 


3r i5- 



B 



12 



..14 



271 



13! ■ 2 ■ C (13) 4tt 2 
(2tt) 13 



13! 



x (2^) 13 12xi 2 



dx 



Hence we finally have 



... 10 691 
i?4 W>TT7TI + 



1051 49 



10 



1051 2494 



llr i3 ' 195r i5 10 0r 14 3r 15 llr 13 100r 14 195r 15 ' 
Plugging this into the right-hand side of (I3.9p . we find that 

23171 1051 16728577 1051 11997107 



T /w 104 164 
T(r) > -r^r^ + 



Let 



then 



1575r 8 1575r 10 80850r 12 200r 13 1576575r 14 800r 15 4204200r 16 ' 
„. , 104 o 164 , 23171 4 1051 , 16728577 , 1051 11997107 



1575 



1575 



80850 



200 



1576575 



800 4204200 ' 



104 

F(r + 2 - 57) -1575 r " 



s 53456 7 12123418 6 1584201713 5 



39375 



984375 



24609375 



32175150922307 4 8482041695506649 3 62685536289049749 2 

151593750000 T + 18949218750000 T + 111718750000000 T 
263340412300075879103 5236655690345652768413 



821132812500000000 



1970718750000000000000 
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Thus T (r) > if r > 2.57. Hence we showed that for r > 2.57, a (r) is a strictly decreasing 
function of r, satisfies the inequality a (r) < g^rzg and lim r ^ +00 a (r) = Numerical eval- 
uation shows that a (2.57) = 0.4709258826 . . . > §§ = 0.4333333333 . . .. From these it follows 
that gjj < a (r) and we conclude the statement of the lemma. □ 

Proof of Lemma \2.2\ To prove the lemma we show that T (r) defined by (|3.7[) is positive if 
< r < 0.3. It is well known (see for example [2]) that 

2n i 1 \« <-//•* + n - ! N 



E7^ ¥ = EH)" 2 I C(2M + 2n-l)r-, , > 0, < r < 1. 

n>l V ' n>0 V 7 

Here C, denotes Riemann's Zeta function. And thus 

S W = E , 2 2 " 2 ,2 = E (-!)" ( 2 " + 2 ) C (2« + 3) r 2n , 

n>l \ n + r ) „> 



E TT^I = E ( n + !)(» + 2) C (2n + 5) r 2 ". 

n>l V ' n>0 

These series are alternating and for fixed < r < 0.3, the sequences 

s n = (2n + 2)C(2n + 3)r 2 ™, 

t n = (n + 1) (n + 2) C (2n + 5) r 2n 
are monotonically decreasing. By Leibniz's theorem we find 

E } 2 2n 2 ^2 = 2C (3) - 4C (5) r 2 + 6C (7) r 4 - 8 r ■ 8( (9) r 6 , 
„>i ( n + r ) 

E TTTT = 2C (5) - 6C (7) r 2 + 12< (9) r 4 - p r • 20( (11) r 6 , 

„>i ( n + r ) 

where < r , ,o r < 1 and < r < 0.3. By plugging these into (|3.7I) . we obtain 

T (r) > (V + 3r 2 + ^ (2( (3) - 4( (5) r 2 + 6( (7) r 4 - 8C (9) r 6 ) 2 

- (6r 2 + 2) (2C (3) - 4C (5) r 2 + 6C (7) r 4 ) 

- (V 2 + (2C (5) - 6C (7) r 2 + 12( (9) r 4 ) + 3 = Q (r) . 

The right-hand side is a polynomial in degree 16. It can be shown that the polynomial 
(r + "y) 16 Q ( f^Joj$ J nas oruv positive coefficients, thus Q (r) is positive on the range < r < 

0.3, hence a is strictly decreasing there. This means that || < a (0.3) < a (r) < a (0 + ) = 6 ^^_ 6 , 
which completes the proof of the lemma. □ 

Proof of Lemma \2.3[ As an application of the Euler-Maclaurin summation formula, Lampret [5] 
proved that 

_ . . v-^ 2n 1 m 3m z — r z 

S (r) = > o H ~ ^ H 5 H 5-+P( TO ; r )) 

^(n 2 + r 2 ) 2 m 2 + r 2 ( m 2 + r 2) 2 6 ( m 2 +r 2)3 

where 

5m 4 + 15m 2 r 2 + 6r 4 

|p(m,r)| < g . 

16 (m 2 + r 2 ) 
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With m = 5we obtain 

_ , . 2n 1 5 75 -r 2 3125 + 375r 2 + 6r 4 ^ , , 

S r > > 7 + — ^ H =■ + t F = Zi r , 

(« 2 + r 2 ) 2 25 + r 2 (25 + r 2 ) 2 6 (25 + r 2 ) 3 16(25 + r 2 ) 5 

and 

Mr)-;^>0,r>0. 

According to ([S3), (ESJ) and ((33]) we have 

(l2r 4 + 2r 2 - l) Z x (r) - 12r 2 

a (r ) < — — — ■ r > 

[ ) 12- (12r 2 + 6)Zi(r) ' 

The right-hand side is a rational function of r. Let 

d (I2r 4 + 2r 2 - l) Z x (r) - 12r 2 _ Vl (r) 

dr~ 12 - (12r 2 + 6) Z x (r) ~~ Vi (0 ' 
where <y9i (r) and "0i (r) are polynomials in r, ( r ) > 0. ip% (r) has degree 45 and 
(r + 9^f) 45 <Pi ^ r -)-i/2 57 ) < if r > since its every coefficient is negative. Hence 

(I2r 4 + 2r 2 - l) Z x (r) - 12r 2 



CO) 

6C(3)-6' 



12 - (12r 2 + 6) Zi (r) 

is a strictly decreasing function of r if < r < 2.57. Its value at 0.3 is 0.9596637512 . . . < 
thus a (r) < gg7^~g if 0.3 < r < 2.57. If m = 4, Lampret's formula gives 

_ , , 2n 1 4 48 -r 2 1280 + 240r 2 + 6r 4 „. . 

Z 2 r = > ? + — H =• H * H = > Sir). 

^i( n2 +r 2 ) 16 + r (16 + r 2 ) 2 6(16 + r 2 ) 3 16(16 + r 2 ) 5 

Since S (r) > r -i+\/2 1 according to (|3.2[) . (|3.3p and (|3.4p we obtain 

(I2r 4 + 2r 2 - l) Z 2 (r) - 12r 2 

a (r) > ' r > 

[ ) 12- (12r 2 + 6)Z 2 (r) ' 

Again, the right-hand side is a rational function of r, let 

d (I2r 4 + 2r 2 - l) Z 2 (r) - 12r 2 _ tp 2 (r) 
d7 12 - (12r 2 + 6) Z 2 (r) ~ ^2 (r) ' 

where <£> 2 (r) and ^2 ( r ) are polynomials in r, ip2 (f) > 0. <p 2 (r) has only negative coefficients, 
hence 

(I2r 4 + 2r 2 - l) Z 2 (r) - 12r 2 



12 - (12r 2 + 6) Z 2 (r) 

30' 



is a strictly decreasing function of r if r > 0. Its value at 2.57 is 0.4360975104. . . > it, thus 



a ( r ) > M ^ ^-3 < r — 2-57 and this completes the proof. □ 

4. Proof of Theorem 11.11 
From the lemmas we conclude that 

a ( + oo).H <a(r)< _C|__ ( o + ), 
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if r > 0. By the definition of a (r) this is equivalent to 

1 



— <S(r)< 



We have 




thus the lower bound is sharp. From Lemma 



CM 1 1 ' ' ' 



r 2 6r 4 30r 6 42r 8 30r 10 66r 12 
asr-> +oo. Plugging this into (I3.ip yields 

, . 13 104 592 404032 

a ( r ) 1 1 1 1 

30 525r 2 2625r 4 1010625?' 6 
as r — > +oo. Hence we obtain a new asymptotic expansion to Mathieu's series: 

S (r) 

1 4r 2 + 1 I 13 104 592 404032 

r 2 + 



2 12 I 30 525r 2 2625r 4 1010625r 6 J 
holds as r — > +oo. This shows that the upper bound is sharp too. 
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